Experimental studies of QCD using flavour tagged jets from DELPHI 
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Identified bbg and qq~/ events from DELPHI are used to measure the ratio of the mean charged particle multiplicity 
distribution between gluon and quark jets. The dependence of this ratio with the jet energy is established using 
about three million Z° decays. Results from all other detectors are discussed and compared. A nice agreement is 
found among all them. The ratio between the normalized total three-jet cross sections of bbg and qqg, q = u,d,s 
events is also determined. The preliminary value obtained indicates that 6 quarks are experimentaly seen to 
radiate less than light quarks due to their higher mass. The suggested experimental error is ~300 MeV for the b 
mass determination at the Mz scale. 



1. INTRODUCTION 

In Quantum Chromodynamics (QCD), quarks 
(q) and gluons (g) are coloured objects that carry 
different colour charges. Quarks have a sin- 
gle colour index while gluons are tensor objects 
carrying two colour indices. Due to this fact, 
quarks and gluons differ in their relative coupling 
strength to emit additional gluons, and, in con- 
sequence, jets originating from the fragmentation 
of energetic quarks and gluons are expected to 
show differences in their final particle multiplici- 
ties, energies and angular distributions. 

The masses of quarks are also fundamental pa- 
rameters of the QCD lagrangian not predicted by 
the theory. The definition of the quark masses 
is however not unique because quarks are not 
free particles and various scenarios are possible. 
The perturbative pole mass, M 9 , and the run- 
ning mass, nig, of the MS scheme are among the 
most currently used. At first order in a s the pre- 
dicted expression for an observable is not able to 
resolve the mass ambiguity and, only when sec- 
ond or higher order terms are included, the mass 
definition becomes known. At orders higher than 
one the renormalization scheme used as the base- 
line of the calculation has to be chosen and this 
contains the information about the mass defini- 
tion. Earlier calculations of the three-jet cross 
section in e + e~ including mass terms already ex- 



ist at 0(a s ) jl],H and have been used to evalu- 
ate mass effects for the 6-quark when testing the 
universality of the strong coupling constant, a s . 
They could not however be used to evaluate the 
mass of the 6-quark, m/,, because these calcula- 
tions are ambiguous in this parameter. Recently, 
expressions at O(af), for the multi-jet produc- 
tion rate in e + e~ are available H and, thus, they 
enable measuring nif, in case the flavour indepen- 
dence of a s is assumed and enough experimental 
precision is achieved. 

There are well known existing difficulties to 
measure all the above parameters in quantitative 
agreement with the predictions from perturbative 
QCD, since partons, quarks and gluons, are not 
directly observed in nature and only the stable 
particles, produced after the fragmentation pro- 
cess, are experimentally detected. However, the 
massive statistics and improved jet tagging tech- 
niques available at LEP presently allow overcomig 
these difficulties by applying restrictive selection 
criteria which lead to quark and gluon jet sam- 
ples with high purities. The selected data samples 
are almost background free and small corrections 
to account for impurities are needed. A smaller 
model dependence than ever is now achieved, 
bringing the possibility to perform quantitative 
studies of quark and gluon fragmentation accord- 
ing to perturbative QCD. 

The analyses reported in here include more 



than 3 million Z° decays as collected by DEL- 
PHI at center-of-mass energies of y/s «Mz. In 
the first analysis, the ratio between the gluon 
jet multiplicity and the quark jet multiplicity, 
r = (N g )/(N q ), is presented and discussed in 
comparison with other detector results. In the 
second study, preliminary values of errors associ- 
ated to the determination of mf, at the Mz scale 
are given. 

2. EVENT SELECTION 

Gluon and quark jets were selected using 
hadronic three-jet events. Jets were mainly recon- 
structed using the Durham algorithm although 
the Jade algorithm was also used in particu- 
lar, to observe the effects due to different angular 
particle acceptance of the various algorithms. 

In the gluon splitting process (g — > qq), the 
heavy quark production is strongly suppressed 
||. Gluon jets can thus be extracted from qqg 
events by applying b tagging techniques. The 
two jets which satisfy the experimental signa- 
tures of being initiated by b quarks are associated 
to the quark jets and the remaining one is, by 
definition, assigned to be the gluon jet without 
any further requirement. Algorithms for tagging 
b jets exploit the fact that the decay products of 
long lived B hadrons have large impact param- 
eters and/or contain inclusive high momentum 
leptons coming from the semileptonic decays of 
the B hadrons. Gluon purities of 94% and 85% 
are achieved when using these techniques, respec- 
tively. Obviously, the quark jets belonging to 
these events cannot be used to represent an un- 
biased quark sample. Thus the quark jets whose 
properties are to be compared with the gluon jets 
must be selected from other sources which in any 
case should preserve the same kinematics. Two 
possibilities have been proposed in the current 
literature. One consists in selecting symmetric 
three-jet event configurations [|||J^1 in which one 
(Y) or the two (Mercedes) quark jets have sim- 
ilar energy to that of the gluon jet. The quark 
jet purities reached are ^52% and ~66%, for Y 
and for Mercedes events, respectively. In a second 
solution {IHH radiative qqj events are selected, 
allowing a sample of quark jets with variable en- 



Table 1 



The three-jet event samples and their correspond- 
ing energy intervals as used in the analysis. 



Event type 


# events 


Jet ener 


gy range 


<W 


2,237 


7.5 GeV - 


42.5 


GeV 


{uds)(uds)g 


552, 645 


7.5 GeV - 


42.5 


GeV 


bbg 


104, 081 


7.5 GeV - 


42.5 


GeV 


Y 


74, 164 


19.6 GeV 


- 28.8 


GeV 


Mercedes 


9,264 


27.4 GeV 


- 33.4 


GeV 



ergy to be collected. In this latter case, misiden- 
tifications of 7's due to the n° background and 
radiative t + t~7 contamination give rise to quark 
jet purities of ^92%. This method gives a higher 
purity but unfortunately suffers from the lack of 
statistics. 

The 6-quark purity in the bbg sample reached 
in the DELPHI analyses is ~93% and for the light 
wds-quarks is ~80%. Table | summarizes the 
number of events selected and their correspond- 
ing energy intervals. 

3. MULTIPLICITIES OF QUARK AND 
GLUON JETS 

Results on the charged multiplicity of quark 
and gluon jets ||,f7| using symmetric Y configu- 
rations and reconstructed with Durham at 24 
GeV gluon jet energy, give a ratio of r k 1.23 ± 
0.04(stat.+syst.) which does not depend on the 
cut-off parameter (y cu t) selected to reconstruct 
jets ||. It is significantly higher than one, which 
indicates that quark and gluons in fact fragment 
differently, but it remains far from the asymp- 
totic lowest order expectation of Cf/Ca = 9/4, 
suggesting that higher order corrections and non- 
perturbative effects are very important to under- 
stand the measured value. A next-to-leading or- 
der correction || in MLLA (Modified Leading 
Log Approximation) at O(y / o7) already lowers 
the prediction towards r values slightly below two 
and exhibits a small energy dependence due to 
the running of a s . However this is still insuffi- 
cient to explain the value of r determined by the 
experiments. Solutions based on the Monte Carlo 



method give a better approximation y . The par- 
ton shower option of the Jetset generator |[o) 
which uses the Altarelli-Parisi splitting functions 
for the evolution of the parton shower reduces 
the theoretical prediction Q for r. At parton 
level, at 24 GeV jet energy, the expected value 
is ~1.4 and it is further reduced to ~1.3 if the 
value of r is computed after the fragmentation 
process. In both cases there is a clear depen- 
dence of r with the jet energy [Q which can be 
parametrized using straight lines with slopes of 
Ar/AE = (+90±3(stat.))-10" 4 GeV" 1 at parton 
level and Ar/AE = (+76±2(stat.)) • 1CT 4 GeV -1 
after fragmentation. The absolute value of r 
predicted at parton level is however largely af- 
fected by the choice on the Qo parameter (cut- 
off at which the parton evolution stops) but has 
negligible influence on its relative variation with 
the energy, i.e., the slope. The DELPHI analy- 
sis uses symmetric and non-symmetric three-jet 
event configurations with quark and gluon jets 
of variable energy, allowing thus all these prop- 
erties and predictions to be tested. A value of 
r = 1.23 ± 0.03 (stat.+syst.) is measured corre- 
sponding to an average jet energy of ~27 GeV. 
The energy dependence of r is also suggested 
at 4<t significance level, with a fitted slope of 
Ar/AE = (+104±25(stat.+syst.))-10- 4 GeV -1 . 

In a recent review fllf all published data from 
various experiments ||jg]ff jlj|[l4fl were used to 
perform a general study of r as a function of 
the jet energy. At present, more data can be 
added to this comparison. These are the new 
analysed DELPHI data sample presented above 
and the most recent measurements of r per- 
formed by CLEO H and OPAL |§ at 4-7 
GeV and 39 GeV average jet energies, respec- 
tively. The updated new DELPHI analysis in- 
corporates two times more statistics than the 
previous analysis Q, therefore, significantly re- 
duces the statistical errors. The analysis from 
CLEO compares the charged particle multiplic- 
ity in T(15) — > ggj decays to that observed 
in e + e~ — * qq~/ just in the continuum. This 
study does not rely on the Monte Carlo simu- 
lation to associate the final hadrons to the ini- 
tial partons and can consequently be fairly con- 
sidered as being model independent. The ob- 
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Figure 1. (a) Mean charged multiplicity of quark 
and gluon jets and (b) multiplicity ratio r as a 
function of the jet energy 



is r = 1.04 ± 0.05 (stat.+syst 



tained value 

The OPAL analysis uses a new technique 
which selects gluon jets at ~39 GeV by divid- 
ing the events into two hemispheres. While one 
of these hemispheres is required to contain two 
tagged quark jets, the other is left untouched be- 
ing regarded as the gluon jet. The result from 
OPAL, expressed for only light uds-quarks, is 
r u da = 1-55 ±0.07 (stat.+syst.). As it can be ob- 
served in figure [|.b all these data agree with the 
predicted energy behaviour of j|,[l2]|ll]] when the 
correction to the quark multiplicity to account for 



the same flavour composition is applied. In our 
case it is: 56% uds's, 33% c's and 11% 6's. The 
OPAL number considering this quark mixture be- 
comes r = 1.48 ± 0.07 (stat.+syst.). 

All these results thus give evidence for an en- 
ergy dependence of r. The measured increase is 

= (+110 ± 13 (stat.+syst.)) • 10~ 4 GeV" 1 , 

representing a ~8ct effect. 

The measured value of r remains systematically 
lower than the Jetset prediction over the whole 
energy range, having an average value of 

r = 1.23 ± 0.01 (stat.) ± 0.03 (syst.), 

which corresponds to an average energy of ^23 
GeV. This ratio can be further expressed as 

r uds = 1.30 ± 0.01 {stat.) ± 0.04 {syst.), 

if r is computed only for the light uiis-quarks, 
extracting the b and c quark contribution to the 
quark jet multiplicity. 

The absolute value of r depends on the re- 
construction jet algorithm. For both the Jade 
and Cone schemes different results are obtained 
w.r.t. the Durham scheme ||||. This is due 
to the combined effect of the different sensitiv- 
ity of the various jet reconstruction algorithms 
to soft particles at large angles and of the ex- 
pected different angular and energy spectra of 
the emitted soft gluons in the quark and gluon 
jets. A precise deconvolution of both effects is, 
at present, impossible iQ. This jet algorithm 
dependence of r becomes however less apparent 
as the jet energy increases. The results from 
OPAL @, r = 1.48 ± 0.07 (stat.+syst.) and 
those from DELPHI |j at ~40 GeV presented 
in this conference, r — 1.43 ± 0.07 (stat.+syst.) 
for Durham and r = 1.52 + 0.11 (stat.+syst.) for 
Jade, agree within errors for the various methods 
and algorithms used. For the low energy inter- 
val, the Jade and Durham jet algorithms give 
a different description of the gluon jet properties 
Q , although the Durham algorithm is in better 
agreement to those, model independent, results 
obtained by CLEO. Hence, the Durham jet al- 
gorithm seems to be better suited to decribe the 



intermediate energy region than the Jade algo- 
rithm is. 

The interpretation of these results in combi- 
nation with those obtained by OPAL jllj and 
ALEPH [^0| restrict the validity of the statement 
that gluon and 6-quark jets have similar proper- 
ties to the jet energy interval around 24 GeV and 
cannot be applied to the whole jet energy spec- 
trum. 

4. GLUON RADIATION IN 6-QUARKS 

For many observable quantities at LEP ener- 
gies, y^s^Mz, quark mass effects usually appear 
in terms proportional to nr^/M^. This represents 
a ^3% correction for m 9 =m(, which in most of 
the cases can be savely neglected. This argument, 
for instance is true for the total hadronic cross 
section Q] but cannot be applied for the differ- 
ential multi-jet cross sections that depend on the 
jet-resolution parameter, y c . The reason being 
the new scale, E c — Mz-J%k, introduced in the 
analysis by the new variable which enhances the 
mass effects in the form m\jE^. = {ml/M^)/y c . 
At yfs kM.z the three-jet production rate for b- 
quarks is in fact suppressed by a factor ~ 5 — 10% 
w.r.t that of light quarks [@,|]|2l]] . This difference 
can then be expressed as a function of mj, [|j and, 
therefore, used to measure its value. 

The experimental observation of such effects is 
however difficult and delicate because the effect 
is after all small and furthermore the correct the- 
oretical framework to resolve the mass definition 
ambiguities is needed. This means that the ob- 
servable has to be calculated including mass ef- 
fects at 0{al). For this purpose a recent cal- 
culation H of the ratio of the normalized three- 
jet cross sections between &-quarks and light uds- 
quarks 
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bd 
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has been performed 

The normalization in R^ d to the total decay 
rates is introduced to cancel possible weak cor- 
rections depending on the top quark mass p2| 
and the ratio of the three-jet cross sections be- 
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Figure 2. Relative systematics uncertainties in Figure 3. R^ d distribution 

R^ d due to fragmentation 



tween b and light uds-quarks minimizes uncer- 
tainties due to the hadronization process. In fig- 
ure ^| the dependence of these uncertainties w.r.t 
the y cut is shown and seen not to exceed 3%o for 
large enough values of ycut- The R^ d distribu- 
tion corrected for detector and fragmentation ef- 
fects is also displayed in figure [| The solid curves 
drawn in the figure are the theoretical 0(a s ) pre- 
diction in steps of 1 GeV. The values of m;, used 
to produce these curves are meaningless since 
they correspond to a calculation at 0(a s ). They 
can nevertheless be used to evaluate the experi- 
mental precision assuming the difference between 
the theoretical curves remains similar to that at 
0(o!g). As can be observed the experimental er- 
ror corresponds then to approximately 300 MeV 
for reasonably high values of y cu t- 
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